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We investigate the scattering off dilute magnetic impurities placed on the surface of three- 
dimensional topological insulators. In the low-temperature limit, the impurity moments are Kondo- 
screened by the surface-state electrons, despite their exotic locking of spin and momentum. We 
determine signatures of the Kondo effect appearing in quasiparticle interference (QPI) patterns as 
recorded by scanning tunneling spectroscopy, taking into account the full energy dependence of the 
T matrix as well as the hexagonal warping of the surface Dirac cones. We identify a universal energy 
dependence of the QPI signal at low scanning energies as the fingerprint of Kondo physics, markedly 
different from the signal due to non-magnetic or static magnetic impurities. Finally, we discuss our 
results in the context of recent experimental data. 
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I. INTRODUCTION 

Topological insulators (TIs) in both two and three 
spatial dimensions constitute an active topic of current 
condensed matter research. 1 7 The non-trivial bulk band 
topology of three-dimensional (3D) strong topological in- 
sulators causes the crossing of surface states at time- 
reversal invariant points in the surface Brillouin zone 
and gives rise to a two-dimensional surface metal. In 
the vicinity of such crossing points, the effective surface 
theory takes the form of a Dirac equation of massless 
fermions, where spin and momentum are locked together. 

A fundamental property of this "helical" surface metal 
is a suppression of backscattering: Electrons with op- 
posite momenta have orthogonal spin projections, such 
that impurity scattering k <o> — k is impossible with- 
out a spin flip. As a result, the metallic state is pro- 
tected from the influence of non-magnetic disorder, and 
weak localization is replaced by weak antilocalizationP^ 
This scenario of forbidden backscattering has been tested 
in recent experiments^"^ utilizing powerf ul Fo urier- 
transform scanning tunneling spectroscop y 16 * 17 ^ (FT- 
STS). In this technique, energy-dependent spatial vari- 
ations of the local density of states (LDOS) are analyzed 
in terms of quasiparticle interference (QPI), i.e., quasi- 
particle scattering processes due to impurities. The QPI 
results obtained on 3d TIs such as Bii-^Sb^ and Bi2Te3 
were found to be consistent with a heuristic picture of 
electron scattering in a helical liquid, with backscatter- 
ing being suppressed. 

These results prompt the question as to how scattering 
from magnetic impurities on the surface of TIs is mani- 
fest in observables such as the QPI patterns obtained by 
FTSTS. In fact, recent experiments 13 on Bi2Te3 doped 
with dilute magnetic Fe atoms purport to demonstrate 
from the QPI pattern signatures of time-reversal symme- 



try breaking. However, one must be careful to distinguish 
a fluctuating magnetic moment from one which is static 
on the large timescale of the STS experiment. The latter 
situation implies magnetic long-range order, whose exis- 
tence requires a sufficient density of magnetic moments 
and low temperature. Then, every impurity moment is 
polarized, and time-reversal symmetry is broken. 18 Inter- 
estingly, it has been shown theoretically that such static 
magnetic impurities do not lead to backscattering being 
visible in QPI; within lowest-order Born app roxim ation, 
a static local field is entirely invisible in QPlJ!212Ql 

In this paper, we focus instead on the case of fluctu- 
ating magnetic impurities, relevant to the dilute limit. 
The interaction between the impurity moment and the 
electrons of the surface metal leads to mutual spin flips, 
such that backscattering could be allowed although time- 
reversal symmetry remains unbroken. It is such spin- flip 
processes which lead to Kondo screening of the impurity 
moment in standard metalsP^ Therefore the key ques- 
tion, also relevant to the experiments of Ref. [I3j pertains 
to the signatures in QPI of the Kondo interaction be- 
tween the helical metal and the impurity. To answer 
this, we solve the problem of a single Kondo impurity on 
the surface of a 3D TI numerically exactly, and calculate 
the induced QPI pattern which of course now includes 
inelastic scattering off the magnetic moment. 

Our main findings are as follows, (i) The magnetic im- 
purity is described by a standard SU(2)-symmetric impu- 
rity model, despite spin-momentum locking and hexago- 
nal warping effects of the surface states of a real TI. As a 
result, the impurity moment is always Kondo screened in 
the low-temperature limit, unless the chemical potential 
is tuned exactly to the Dirac point, (ii) While scattering 
off a Kondo impurity does not open new scattering chan- 
nels in momentum space as compared to a non-magnetic 
impurity, it leads to a distinct energy dependence of the 
QPI pattern, which moreover exhibits universal scaling 
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in terms of both scanning energy and temperature. A 
strong enhancement of the QPI intensity near the Fermi 
level is therefore a signature of scattering caused by fluc- 
tuating magnetic impurities. 

The body of the paper is organized as follows. We 
start by introducing the model and methods in Sec. [TTJ 
The Kondo effect on the surface of 3D TIs is discussed in 



(a) 



Sec. Ill Sec. IV is then devoted to the QPI patterns from 
Kondo impurities, with an emphasis placed on universal 
features. The QPI signal from non-magnetic impurities 
is shown for comparison in Appendix [B] while we briefly 
discuss a static magnetic impurity in Appendix [C] The 
implications of our results for experimental data are dis- 
cussed in the concluding section [V) 

We note that the Kondo effect on the surface of 3D 
TIs was discussed before in Refs. 22 22 



but without 

taking into account hexagonal warping and without a 
discussion of QPI. QPI patterns for surfaces of 3D TIs 
have been calculated for different types of impurities in 
Refs. 110113)19)20(241 but no link to Kondo physics was 
made. Very recently, inelastic scattering from excited 
states of magnetic impurities was discussed in Ref. [25j 
leading to features at elevated-energy in tunneling spec- 
tra. 



II. MODEL AND METHODS 
A. Effective surface metal 

Surface states of 3D TIs are described by an effective 
Dirac theory. However, such a linearized model applies 
only in the immediate vicinity of the crossing point of 
the surface bands, while lattice effects must be taken 
into account at higher energies. For Bi2Te3 this re- 
sults in a breaking of the continuous rotation symmetry 
around the Dirac point down to C$ v , leading to so-called 
hexagonal war ping of the iso-energy contours,^ as is seen 
experimentallyP^ 3 | 27 | 2S | 

The free Hamiltonian of the surface metal reads^^U 



H = J d 2 k 



kf 14 



kf 



(i) 



where 



% k = hv F [(k x a) • e z + A 2 k 3 cos(3</> k )a 2 ] - /i . (2) 

Here a is a vector of the Pauli matrices, k = |k| is the 
magnitude of the momentum vector relative to the Dirac 
point, and (/>k = tan -1 (k y / f k x ) is its azimuthal angle mea- 
sured with respect to the x axis. The Y—K direction 
thus corresponds to 0k = while Y-M corresponds to 
= tt/6, following Ref. [24]. The cubic term ex A 2 ac- 
counts for hexagonal warping, and \i denotes the chemical 
potential. 

The spectrum of the above Hamiltonian (with h = 1 
hereafter) is given by 



E±(k) = ±v F \Jk 2 + [A 2 k 3 cos (30 k )] z - p 



(3) 
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FIG. 1: (color online) (a) Density of states, p^°\uj), for the 
effective surface theory, Eq. ^ in the absence of impurities, 
for A 2 = 2.23. At energies around uj — —p one obtains linear 
pseudogap behavior, which is modified at higher energies due 
to hexagonal warping, (b) QPI pattern Ap(q, uj) for a TI with 
p(°\cj) as in (a), using T°(cj) — 1 appropriate for a potential 
scatterer in Born approximation. Plotted for /jl/vf = 0.137 = 
100 meV and scanning energy uj — 200 meV. 



free Green function, G^°)(k, uj) = [uj + i0 + — 
1 , takes a diagonal form in the quasiparticle basis, 



The 
#k]" 

G^(Kuj) = S ab /[oj + iO+-E a (k)], a,6 : 
of states (DOS) follows from 



d=. The density 



~ImTte (0) (k,o;), 



(4) 



where the trace accounts for the sum over ± as well as 
k, and TV is a suitable normalization factor (equal to 
the number of k points).^ p(°\uj) is linear in |cj| at low 
energies around uj = — /i, characteristic of massless Dirac 
fermions. 

In the following we shall employ parameters vp/a^ = 
0.73 eV (ao = 1 is a lattice constant acting) as our energy 
unit, and A 2 = 2.23 to make contact with experiments.^ 
For convenience we shall use vp/ao as high-energy cut-off 



for the conduction band, p(°\uj) = p^°\uj)0(vf - 



-Ml) 



this is mainly needed to generate an input for the nu- 



merical treatment of the impurity problem in Sec. Ill 



(For a 3D TI, a natural cutoff is set by the size of the 
bulk gap.) Generically, the Fermi level is not at the Dirac 
point; for example, in Ref. [13] the chemical potential is 
p/vp = 0.137 = 100 meV. Below we consider this case 
explicitly, and also the special case where p = 0, which 
is potentially attainable, since surfaces of TIs can be in- 
dividually gated. 

The local Green function on the surface of the TI has 
a matrix structure in spin space which turns out to be 
diagonal: 



G<°> (r = 0,o; n ) =!/„„,„, 
in terms of Matsubara frequency uo n . Here 

I „wwy, , ' P 

Jc 



UJ n ,/J, 



I 



Vpk 2 



vjA^k 6 cos (30 k )' 



(5) 



, (6) 



where / k = / j^s, such that P (0) (lu) = -2Im/ W)J[i /7r. 
As shown in Appendix |Aj the off-diagonal elements of 
G(°)(r = 0) vanish on angular integration. 
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B. Impurities and QPI 

The FTSTS technique exploits scattering of charge 
carriers from imp uritie s in an otherwise translationally- 
invariant system To this end, a real-space map 
of the tunneling conductance is recorded at fixed en- 
ergy. The Fourier transform of this map to momen- 
tum space reveals characteristic wavevectors of LDOS 
inhomogeneities which can be understood as energy- 
dependent Friedel oscillations (or QPI). In the simplest 
approximation, these wavevectors correspond to scatter- 
ing processes of quasiparticles between different points 
of the dispersion iso-energy contour at the scanning 
energy 

If scattering centers are dilute, it is sufficient to con- 
sider a single impurity. Its effect is described by the T 
matrix, Tk k'(co>), such that the full electronic Green func- 
tion reads 

G(k,k» = G(°)(k ) w)<5 k , k / 

+ G (0) (k, w)T k>k , (w)G (0) (k', u) , (7) 

which remains energy-diagonal in equilibrium/linear re- 
sponse. The real-space LDOS is 

p(r,cj) = --Im /e iq r Tr G(k,k-q,cj), (8) 



where the trace again accounts for the sum over k as well 
as the spin components. For isotropic scattering and an 
inversion-symmetric host, the Fourier transform p(q, oS) 
is real, and its impurity-induced piece is related to the T 
matrix via 

Ap(q, u) = - -Im Tr (k, cj)T k>k _ q (cj)G* > (k - q, u>) . 

7T 

(9) 

// the T matrix is diagonal in spin space (see below), 
decomposition using the Pauli matrices gives 



r k>k ' M = K k , (cj)I + T" (uj)cj z 



(10) 



with T^Uuo) = [T k>k /, t (o;) ±T k>k /^(o;)] /2. For TI sur- 
faces in the absence of hexagonal warping it was shown 
in Ref. |20] that only the part proportional to I leads to 
a modulation of the spin-integrated LDOS in Eqs. (8|9). 
(The part proportional to a z causes opposite modulations 
for both spin directions.) Using symmetry properties, we 
have verified that this still holds for the full model includ- 
ing hexagonal warping: the argument parallels that given 
in Appendix [Aj 

Specializing further to a point-like impurity with 
^k,k'M = f(u>)/N, we have 



Ap(q,cj) = - -Im T°(uj) Tr G^(k, ou)G^ \k - q, u) 

7T 

(ii) 

such that the momentum dependence of the QPI signal at 
fixed uj is completely determined by . A sample QPI 



image is displayed in Fig. [I]- this is similar to published 
results. 19 One clearly observes a breaking of the circular 
symmetry due to hexagonal warping of the iso-energy 
contour at this energy. 

As we show below, a (local) Kondo impurity does not 
modify the overall momentum dependence of the QPI 
pattern, but will lead to a non-trivial energy dependence. 



C. Anderson impurity model 

To describe a dynamic magnetic scatterer, we consider 
Anderson's model for a point-like correlated impurity,^ 
H = H + H imp + # hy b, with 

#imp = e d {f^ + hf) + Unfaf , 

H hyh = gJ24**(r = o) + n.c.. (12) 



Here * ff (r = 0) = E k W^, 



and e d < 



and U > —€d are the local level energy and Coulomb 
repulsion, respecively. 

Considering the spin-momentum locking and the 
hexagonal warping of the TI surface electrons in i^o, one 
might have expected an unconventional impurity prob- 
lem. However, a standard SU(2) spin-symmetric impu- 
rity problem is obtained, with the complexity of the heli- 
cal surface metal entering only through the unusual DOS, 
Eq. Q. The derivation of such a pseudogap Anderson 
(or Kondo) model has been established before for im- 
purities in d-wave superconduc tors^ as well as for TIs 
with a perfect Dirac structure, 1 22 1 23 1 and we give here 
an efficient proof which also covers hexagonal warping. 
Rather than using a decomposition into angular modes 
as in Refs. I22|23| a more direct way is to use the path- 
integral formulation. Since the conduction-electron bath 
is Gaussian, we integrate it out exactly^ to derive a lo- 
cal retarded impurity problem. The local non-interacting 
part of the action for the d-levels in Matsubara formalism 
reads 

So = -T^d T [g 2 G^ (r = 0,iu> n ) - I(e d - iu n j\ d , 

(13) 

with G^ (r = 0,iLJ n ) given in Eq. (51) and d = 
The resulting local model is equivalent to the standard 
Anderson model because (r = 0, iuj n ) is diagonal. 
The hybridization function characterizing this impurity 
problem is 



AM 



-Im g 



(14) 



Importantly, the impurity Green function G d (u) is di- 
agonal due to the absence of off-diagonal terms in the 
quadratic action of the d-levels, Eq. (13). In the absence 
of a magnetic field, we thus have 



G d {w) 



1 



w + iO+-e d - 5 2 /^-S d H 



I. (15) 
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Here £ d (co>) is the interaction part of the impurity self- 
energy. The corresponding T matrix is then T^'(uo) = 
g 2 G d (uj)/N oc I, and one expects a non-trivial re- 
sponse from Anderson-type impurities in QPI according 
to Eq. <[TTJ. 



III. KONDO EFFECT 



The Anderson impurity model, Eq. (12), can describe 



both the formation and the subsequent screening of lo- 
cal magnetic moments. Charge fluctuations induce two 
broadened peaks ("Hubbard satellites") in the impurity 
spectral function at and (U + e^). The physics at 
energies smaller than min(— e^, U + e^) is dominated by 
spin-flip processes leading to Kondo screening below a 
temperature Tk - a phenomenon which is sensitive to 
the conduction-band DOS near the Fermi levelP^ Given 
the unusual DOS of the TPs surface metal, which van- 
ishes linearly at uj = —fi, the Kondo effect should be 
discussed separately in the two cases: (A) the chemical 
potential is tuned to the Dirac point, \i = 0; (B) finite 
chemical potential, /i ^ 0. 

We will obtain numerical results for the Anderson 
impurity model using Wilson's numerical renormaliza- 
tion group (NRG) technique.^ In the following, the hy- 
bridization function, Eq. (14), is discretized logarithmi- 
cally using A = 3, and ~ 3000 states are retained at 
each step of the iterative diagonalization. The results of 
z = 3 interleaved calculations^ are th en com bined for op- 
timal results. The full density matrhP^^ is calculated, 
and from it the impurity spectral function is determined 
numerically^ as a full function of energy, cj, at arbitrary 
temperature, T. The real part of G d is obtained using 
Kramers-Kronig relations. The T matrix and hence QPI 
can then be calculated via Eq. (11). 



We note that the full DOS at the surface of a 3D TI will 
also have higher-energy contributions from bulk states, 
not captured by our modelling. Consequently, our cal- 
culation cannot establish a quantitative link between the 
parameters of the Anderson model and Tk. At present, 
there is no experimental information available on the ac- 
tual values of Tk for concrete TI materials and impu- 
rities. Therefore we choose parameters of Eq. (12) such 



that Tk attains values of order 1 K. For most calculations, 
we shall employ parameters U/vp = 0.3, e/vp = —0.1, 
and g/vp — 0.168 (chosen to match some elevated-energy 
features of the data in Ref. 13), corresponding to a mod- 
erately correlated impurity. Results are shown for T = 
unless otherwise noted. 



A. Chemical potential at the Dirac point 



For fi = 0, the density of states at low energies is 
p^(uj) oc \uj\ r with r = 1. The low-energy physics 
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is thus that of the pseudogap Kondo model. 
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FIG. 2: (Color online) Impurity spectral function vfA(uj) 
vs energy uj/vf at T — for U/vf = 0.3, c/vf = —0.1, and 
g/vF = 0.168. (a) Deep in the local moment phase for \i — 0. 
Inset shows a close-up in the vicinity of the Fermi level, with 
the linear soft-gap in the free density of states showing up 
at low energies (red dotted line is vfA(uj) ~ |cj|). (b) Deep 
in the Kondo screened phase, \i — 100 meV. Inset shows a 
close-up of the Kondo resonance. 



the case of particle-hole symmetry in both the impu- 
rity (ed = —U/2) and the bath, screening is absent for 
all parameters, and the impurity moment remains free 
down to the lowest energy /temperature scales. In con- 
trast, if particle-hole symmetry is broken, a quantum 
phase transition occurs, separating the local-moment and 
Kondo screened phases. However, the latter requires 
strong particle-hole asymmetry and impurity-host cou- 
pling, such that screening is less likely to occur at the 
Dirac point. 

A typical impurity spectral function A(uo) = 
— Im G d {uj)/n for /x = in the local- moment phase 
at T = is shown in Fig. [2|a). At high energies 
(uj/vp « —0.1 and uo/vp « 0.2) clear signatures of the 
Hubbard satellites are observed. At low energies (see 
inset) spectral weight is suppressed due to the pseudo- 
gapped free density of statesPHEH 



B. Chemical potential away from the Dirac point 

For finite chemical potential, \i ^ 0, there is a finite 
density of states at the Fermi level. The impurity is 
always screened by the Kondo effect on the lowest en- 
ergy scales (although the Kondo temperature, Tk, itself 
might be very small). Screening is reflected in the im- 
purity spectral function by a narrow resonance around 
the Fermi level. Together with the high-energy Hubbard 
satellites, this three-peak structure is the classic hallmark 
of the Kondo effect. 

In Fig. [2|b) we plot the spectral function for an im- 
purity with for the same parameters as in panel (a), but 
with fi = 100 meV. The inset shows a close-up of the 
Kondo resonance, of width Tk ~ 10~ 4z vp ~ 1 K [we de- 
fine T K via A(uj = T K ) = \A(uj = 0)]. 

For fi values closer to the Dirac point, the impurity 
model will display non-trivial crossover phenomena,^ 
different from those of the standard Kondo problemp^ 
due to the non-constant DOS and the proximity to the 
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/jl = quantum phase transition. Such crossovers can be 
expected when \i < Tk and are not present in Fig. [2jb). 



IV. QPI FROM DYNAMIC MAGNETIC 
IMPURITIES 



We now calculate the QPI pattern, Ap(q, cj), induced 
by a dynamic magnetic impurity on the surface of a 3D 
TI, first for the generic situation of finite chemical po- 
tential and then for the special case where the chemical 
potential is tuned to the Dirac point. We recall that 
Ap(q, uj) is real for our case of a single impurity; experi- 
ments typically extract its absolute value. 



A. Kondo phase 

A finite chemical potential implies Kondo screening at 
lowest temperatures. We divide our analysis into the 
regimes of elevated and low energies. 



200 meV 



220 meV 



240 meV 
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1. Elevated energies, uj ^> Tk 

The QPI pattern obtained at high energies 200- 
360 meV for \i = 100 meV is shown in Fig. [3j Upon 
increasing the scanning energy, the high-intensity peaks 
move outwards and become more pronounced. This is 
to be expected from the increasing diameter of the Fermi 
surface and the increasing importance of hexagonal warp- 
ing, which is due to the underlying lattice structure and 
gives rise to the six-fold symmetry. We recall that our 
modelling neglects bulk bands which will contribute to 
the signal at energies beyond the bulk gap! 13 | 4S | 

Cuts through the QPI along the Y-K and Y-M direc- 
tions shown in Fig. fallow a more detailed analysis of the 
high-energy behavior. One observes that the shape of the 
Hubbard satellites is manifest as a non-montonic inten- 
sity variation of the QPI peaks with energy, in contrast 
to what would be observed for non-magnetic impurities 
(see Appendix|B|. Such a feature is reminiscent of the ex- 
perimental findings in Ref. [I3j We note that, in a more 
complete modelling of the impurity, non-monotonic be- 
havior could also arise from excited crystal-field states of 
the magnetic impurity. 



2. Universal Kondo regime, uj ~ Tk 

The pronounced build-up of impurity spectral intensity 
in a very narrow energy window O(Tk) around the Fermi 
level is the characteristic signature of the Kondo effect 
see Fig.[2jb)]. This results in a similarly characteristic 
evolution of the QPI pattern at low energies. 

For the parameters used in Fig. |3j the Kondo scale is 
Tk ~ 10 -4 vf, and so we probe the system around these 
energies in Fig. [5] The intensity of the QPI peaks indeed 



FIG. 3: (Color online) QPI patterns Ap(q, uj) at different 
energies uj for a dynamic magnetic impurity, with parameters 
A 2 2.23, /x 100 meV, U/v F 0.3, e d /v F -0.1, and 
g/v F = 0.168. 



grows rapidly in this regime due to the Kondo effect. 
This is further highlighted in Fig. [6j where cuts through 
the QPI pattern are shown. (We note the different scales 
in Figs. [3] and [5] as well as Figs. [4] and |6j respectively.) 

At low energies in the Kondo regime, there is a uni- 
versal scaling collapse of the impurity Green function in 
terms of uj/Tk and T /Tk- As in this regime (and in the 
scaling limit) the free Green function (k, uj) is essen- 
tially independent of energy (and strictly independent of 
temperature) one naturally expects universality and scal- 
ing of the entire QPI pattern. As we show in Fig. [7] this 
is indeed the case. 

The T = peak intensity for different choices of the 
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FIG. 4: Cuts through the QPI signal for a dynamic mag- 
netic impurity along (a) the Y-K direction, and (b) the Y-M 
direction, for a system with the same parameters as in Fig. [3] 
Energies shown are from 160 meV-360 meV in steps of 20 meV 
in order of increasing peak position. 
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FIG. 5: (color online) QPI patterns |Ap(q, for a dynamic 
magnetic impurity in the low-energy Kondo regime, with the 
same parameters as Fig. [3] but now with scanning energies 
very close to the Fermi level. Note the different scale as com- 
pared with Fig. [3] also, Ap(q, u) has sign changes, see Fig. [6] 



system parameters, and thus the Kondo temperature, is 
shown in Fig. [7^a). At low energies the data is well de- 
scribed by 
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FIG. 7: (color online) (a) Scaling collapse of the T = QPI 
peak intensity at low energies for different bare parameters 
and hence different Kondo temperatures, Tk. The peak posi- 
tion at q = q p is essentially pinned at low energies (see Fig.|6| 
as /jl/vf = 0.1 is kept fixed. Circles: (U/vf^/vf, 9/vf) = 
(0.2,-0.1,0.126), T K « 10~ 2 K; Squares: (0.2,-0.1,0.168), 
T K « 2K; Diamonds: (0.3,-0.15,0.126), T K « 6 x 
10~ 5 K; Crosses: 0.3,-0.15,0.168), T K « 0.1K. For illus- 
tration we have also plotted (Green Stars) the peak inten- 
sity for (U/vf,€/vf, g/vF) = (0.2,-0.1,0.175) and a differ- 
ent ii/vf — 0.137, yielding Tk ~ 20 K. Red dashed line: 
Low-energy asymptote Eq. |l6| ). Blue dotted line: High- 
energy asymptote Eq. (17). (b) Scaling collapse of the uj — 



peak intensity at low temperatures for the same systems as in 
(a). Red dashed line: Low-temperature asymptote Eq. (18). 
Blue dotted line: High-temperature asymptote Eq. (19). De- 



viations from scaling at higher T become apparent for the 
T K = 20 K data. 



p(q p ,u;,T = 0) 



= 1 - a 



with p(q p ) = p(q p ,w = 0,T = 0) 



, w < T K , 



(16) 



where q p denotes the scattering vector of the peak. The 
quadratic behavior reflects the Fermi liquid nature of the 
Kondo ground state. Similarly, at elevated energies the 
asymptote is given by 



p(q p ,u;,T = 0) 
p(Qp) 



]n(a'\Lj\/T K 



-, T K «w«^, 

(17) 

corresponding to spin-scattering processes in the vicin- 
ity of the local-moment fixed point. These forms reflect 
universality in the scaling limit (Tk —> 0) of the Kondo 
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FIG. 6: Cuts through the QPI signal for a dynamic magnetic 
impurity in the low-energy Kondo regime along (a) the Y-K 
direction, and (b) the Y—M direction, for a system with the 
same parameters as in Fig. 3 Ener gies shown are uj/vf 
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in order of increasing peak intensity. 
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impurity, with universal coefficients a, a' = 0(1). Cor- 
rections arise from two sources: There is a leading linear- 
in-cj term in Eq. (16) because the QPI pattern involves 
both real and imaginary parts of G^(cj), and the real 
part of the self energy is generically non-zero. The linear 
energy variation of contributes to this linear term as 
well (this is where an explicit \i dependence enters), and 
it also causes deviations in Eq. ( fl7| ). These corrections 
are suppressed by the factor cj/min(|/i|, vp) and are small 
for the data shown in Fig. (7^a). 

The temperature scaling of the uo = peak intensity 
is shown in Fig. (7^b). In the low-temperature limit one 
finds 



p(q p ,o; = 0,T) 
P(<Jp) 



T 

5k 



T « T K , (18) 



characteristic of Fermi-liquid behavior. Indeed, we find 
that the Fermi liquid relation a/b = 3/tt 2 is well-satisfied 



from our numerical calculations. Similarly to Eq. ( 17) the 
spin-flip scattering processes in the vicinity of the local- 
moment fixed point lead to the temperature dependence 



p(q p ,uj = 0,T) 



ln(6'T/T K ) 2: 



r K «T«^, (19) 



where 6,6' = 0(1) are again universal. Deviations 
from universality arise as above; they are visible for the 
highest-T K data in Fig.fffb). 

The characteristic feature of the Kondo effect is thus 
the rapid increase of QPI peak intensity as the scanning 
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energy approaches the Fermi level, together with the (ap- 
proximate) scaling collapse in u/Tk and T/T&- These 
features should be readily observable provided that ex- 
periments are performed at temperatures of order Tk or 
lower. 

Finally, we comment briefly on the case where a (small) 
magnetic field acts on the dynamic magnetic impurity. 
Since T^ f (uj) and the impurity model in the absence 
of a field is SU(2) symmetric, it is sufficient to discuss 
the spin-summed impurity spectral function. For small 
fields Q(Tk), this is known to develop a split Kondo 
resonance,^ signatures of which should show up in QPI. 
For larger fields ^> Tk, the Kondo effect is destroyed 
entirely, and only the high-energy Hubbard satellites re- 
main thus leading to a behavior qualitatively similar to 
Fig-i 

B. Local- moment phase 

In the local- moment regime, realized if the chemical 
potential is tuned to the Dirac point /i = (unless the 
Kondo coupling and asymmetry is very strong), there 
is no low-energy Kondo resonance in the impurity spec- 
tral function, although the Hubbard satellites remain [see 
Fig.§a)]. 

As a result, the high-energy features of the QPI pattern 
are essentially the same as those in Figs. [3j [4j but there 
is no build-up of QPI peak intensity in the vicinity of 
the Fermi level. The same applies to the situation with 
non-zero \i but temperatures T ^> Tk- 

V. CONCLUSIONS 

In this work we have studied modulations in the LDOS 
caused by dilute magnetic impurities on the surface of 
3D TIs. Despite the coupling of orbital and spin degrees 
of freedom in the helical surface metal and hexagonal 
warping due to the underlying TI lattice structure, the 
quantum impurity problem itself is of conventional type, 
such that Kondo screening is generically present in the 



low-temperature limit (unless the chemical potential is 
tuned to the Dirac point). 

We identified the energy dependence of the QPI peak 
intensity as the fingerprint of dynamic magnetic impuri- 
ties (as opposed to non-magnetic or spin-polarized impu- 
rites, see appendix). At elevated energies, non-monotonic 
QPI peak intensity for scanning energies is observed due 
to Hubbard satellites in the impurity spectral function. 
At low scanning energies, Kondo screening of the impu- 
rity produces a strong build-up of QPI peak intensity, 
whose energy and temperature dependence are universal 
functions of uo/Tk and T/Tk- 

However, the momentum dependence of the QPI signal 
(at fixed energy) is identical for the different types of 
impurities, due to the fact that only the spin-diagonal 
part of the T matrix produces modulations in the charge 
channel.^ This casts doubts on the interpretation of the 
experimental data given in Ref. [13j where the magnetic 
character of dilute Fe impurities was made responsible for 
the appearance of new QPI wavevectors. At present, the 
source of this QPI signal is unclear, and more systematic 
studies are called for. 
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Appendix A: Local Green function and symmetries 

Here we consider the form of the local free Green func- 
tion, and show explicitly that its off-diagonal components 
vanish. The inverse Green function for Matsubara fre- 
quencies in spin-space reads 



G\iu n ^ = ( l ^ + »-f k \™\^ . - VF i% + J kx) ^ ) (Al) 
v n ' 7 y -v F (k y -ik x ) iuj n + fj, + A z k 6 cos 30 k / 



where the local Green function assumes the form 
G°(iuj nj r = 0) = 



kdk f dcj) 1 / iuj n + ji - A 2 k 3 cos 30 v F (k sin + ik cos 0) 

~27 J 2^ (iu n + /i) 2 - v F (k 2 + A A k 6 cos 2 30) V v f (k sin - ik cos 0) iuo n + /i + A 2 k 3 cos 30 

(A2) 

The off-diagonal terms do not survive the angular integration, such that the remaining diagonal components simplify 
to 

fiOf _ n\ _ tt f f W iUJn +/i / A 

^T-U)-1J ^ J 27r( ^ n + /i) 2_^2 (A ,2 +A 4 A ,6 cos 2 30) ^ 
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a o^r-M 



FIG. 8: Cuts through the QPI signal for a scalar impurity 
along (a) the Y-K direction, and (b) the Y-M direction, for a 
system with fi = 100 me V and a potential scattering strength 
V/vf = 0.1. Energies shown are from 0meV-360meV in 
steps of 20meV in order of increasing peak position (20 meV 
- 0.027vf). 



mixed-valent impurity atom. Then, the impurity spec- 
tral function consists of a single peak centered around 
uj = ed- In contrast to the scalar impurity discussed 
above the T matrix of the resonant level possesses an 
appreciable imaginary part (see Eq. (15) with £ d = 0). 

Cuts through the QPI pattern along the Y-K and T- 
M directions at different scanning energies are plotted 
in Fig. [9] for \i — 100 meV, ed/vp = —0.1 and g/vp = 
0.168. For this resonant level, the intensity of the peak 
is reminiscent of the high-energy response produced by 
a Kondo impurity obtained for U > (simple charge 
fluctuations are responsible for the evolution of the QPI 
signal at high energies in both cases). However, crucially 
there is no build-up of QPI intensity at low energies, since 
there is no Kondo effect (see by contrast Fig. [6|. 



This immediately allows one to read off the energy eigen- 
values in Eq. 

Appendix B: QPI from non- magnetic impurities 

For completeness, we show here the QPI signal of non- 
magnetic impurities, which has largely been calculated 
and discussed in Refs. I10ll3ll9l20l24l 



Appendix C: QPI from static magnetic impurities 

The simplest model for a static, i.e. polarized, mag- 
netic impurity is a local magnetic field, corresponds to 
a spin-dependent version of the potential scattering case 
considered above. Thus we have H = Hq + Hh, with 



1. Scattering potential 



H h = h^^Ar = 0)<x z * ff (r = 0) 



(CI) 



The TI surface metal with a non-magnetic scalar impu- 
rity or potential scatterer is described by H = Ho + Hy , 
with 



H v = vJ2H(r = 0)* <7 (r = 0) , 



(Bl) 



Within the Born approximation, the T matrix is then 



N 



(C2) 



The exact T matrix is 
Tk, k '(w) 1 



V 



N 1 - Vf u , t 



i=-tV)i 



(B2) 



with given in Eq. (|6|, while lowest-order Born ap- 
proximation (valid for small V) corresponds to T°(uj) = 
V which is real and constant. The QPI pattern in this 
approximation is displayed in Fig. [TJb). 

Since impurities are often not weak, the Born approxi- 
mation may be insufficient. Results using the full T ma- 
trix, which now includes a small imaginary part, are plot- 
ted in Fig. [8] for V/vp = 0.1. We observe that the inten- 
sity of the peak along the Y-K direction simply increases 
monotonically as the scanning energy \u\ increases (for 
both uj > and uj < 0). 



2. Resonant level 

A so-called resonant level impurity is obtained by set- 
ting U = in Eq. (12), physically corresponding to a 



i.e. T 1 ^) = h, such that there is no response in QPI 
due to Eq. (11). Beyond the lowest-order Born approxi- 
mation, a (weak) response similar to that of a potential 
scatterer is induced, see Refs. H8ti20l 
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FIG. 9: Cuts through the QPI signal for a resonant level 
along (a) the Y-K direction, and (b) the Y-M direction, for a 
system with \i — 100 meV, U = 0, € c i/vf = —0.1 and q/vf — 
0.168. Energies shown are from 0meV-360meV in steps of 
20 meV in order of increasing peak position. 
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